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Abstract. The second pulsar catalogue of the Fermi Large Area Telescope (LAT) will contain 
in excess of 100 gamma-ray pulsars. The light curves (LCs) of these pulsars exhibit a variety 
of shapes, and also different relative phase lags with respect to their radio pulses, hinting at 
distinct underlying emission properties (e.g., inclination and observer angles) for the individual 
pulsars. Detailed geometric modelling of the radio and gamma-ray LCs may provide constraints 
on the B-field structure and emission geometry. We used different B-field solutions, including 
the static vacuum dipole and the retarded vacuum dipole, in conjunction with an existing 
geometric modelling code, and constructed radiation sky maps and LCs for several different 
pulsar parameters. Standard emission geometries were assumed, namely the two-pole caustic 
(TPC) and outer gap (OG) models. The sky maps and LCs of the various B-£eld and radiation 
model combinations were compared to study their effect on the resulting LCs. As an application, 
we compared our model LCs with Fermi LAT data for the Vela pulsar, and inferred the most 
probable configuration in this case, thereby constraining Vela’s high-altitude magnetic structure 
and system geometry. 


1. Introduction 

Pulsars are considered to be cosmic lighthouses that rotate at tremendous rates and are highly 
magnetized neutron stars (NS) [T]. The fact that pulsars are embedded in such extreme 
conditions make them valuable laboratories for studying a wide range of topics, including: 
nuclear physics, plasma physics, electrodynamics, magnetohydrodynamics (MHD), and general 
relativistic physics [2]. Pulsars emit radiation across the electromagnetic spectrum, including 
radio, optical. X-ray and gamma ( 7 ) rays [1]. We focus on 7 -ray pulsars, specifically the Vela 
pulsar, which was detected [3] by the Fermi Large Area Telescope (LAT) [3]. The Vela pulsar is 
the brightest persistent GeV source. Fermi was launched in 2008 and has discovered in excess 
of 100 7 -ray pulsars. The LAT has a very large field of view of 2.4 sr which enables it to observe 
20 % of the sky at any instant, scanning the entire sky in a time frame of a few hours. 

1.1. Geometric pulsar models 

There are several models which can be used for the modelling of high-energy (HE) emission from 
pulsars. These models include the two-pole caustic (TPC) (the slot gap (SG) [S] model may 



Figure 1. A schematic representation of geometric pulsar models. The TPC emission region 
extends from i?NS (NS radius) up to i?LC (light cylinder radius), OG region from i?NCS (null 
charge surface radius) to i?LC) and the PSPC from i?NS to Ri^c [covering the full open volume 
region]. 


be its physical representation) model [B], outer gap (OG) model [7],[S] and pair-starved polar 
cap (PSPG) model [B]. Figure 1 illustrates these geometric pulsar models, and their emission 
regions [B]. 

Gonsider an (fl, fl) plane, with fl (the magnetic moment) inclined by an angle a with respect 
to the rotation axis n (the angular velocity). The observer’s viewing angle ( is the angle between 
the observer’s line of sight and the rotation axis. The gap region is defined as the region where 
the relativistic particles originate and particle acceleration takes place. The emissivity of HE 
photons within this gap region is assumed to be uniform in the co-rotating frame and the 7 -rays 
are expected to be emitted tangentially to the local magnetic field in this frame m- The gap 
region for the TPG model extends from the surface of the NS along the entire length of the 
last closed magnetic field lines, up to the light cylinder, as indicated by the dashed lines in 
Figured) For the OG model, the gap region extends from the null-charge surface (NGS), where 
the Goldreich-Julian charge density is pcj = 0 m, to the light cylinder, as indicated by the 
shaded region (the emission region may be located at slightly smaller co-latitudes compared to 
the TPC). The PSPC is the gap region that extends from the surface of the NS to the light 
cylinder over the full open volume [9]. 

1.2. Assumed magnetospheric structure 

The magnetospheric structures studied in this paper include the static m and retarded vacuum 
dipole m- The (aligned) static dipole is a special case of the retarded dipole and is described 
by the following B-field equations in terms of spherical coordinates in the laboratory frame: 
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The retarded dipole is described by the following B-field equations M- 
Bret,r = [cos a COS 6 + sin a sin 0(rn sin A + cos A)] 

Bret,4> = “ S™ “[(^n “ 1) COS A] 

Bret,e = (cos a sin 6 + sin a cos 6 [—rn sin A + (r^ — 1) cos A]) 
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By setting equal to zero the retarded field simplifies to the general (non-aligned) static 
dipole, with r the radial distance. The static dipole field is studied for numerous reasons. Two 
of them are: (1) calculations are simpler for this B-field, and (2) when the results for the static 
dipole are compared to those for the other B-fields, the importance of the near-iiLc distortions 
in the B-fields for predicted radiation characteristics can be gauged m- 

In this paper we will study the impact of different magnetospheric structures on the 
predictions of 7 -ray pulsar LCs. The layout is as follows: §2 describes the method we used 
to construct sky maps, LCs, and contour plots for the different combinations of the two 
B-fields and two geometric models. Section 3 contains our results, §4 our discussions and §5 
contains the conclusions and future aims. 


2. Method 

We used an existing geometric modelling code [TO] in which different B-field solutions and 
geometric models are implemented. We constructed sky maps, which are defined as the intensity 
per solid angle as a function of phase and and LCs for the B-field and radiation model 


Table 1. Best-fit (a, C) values for the Vela pulsar. 
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combinations, using a 5° resolution for a and C,. LCs are obtained by making a constant-^^ cut 
through each sky map. 

After the preparation of the sky maps and LCs, a statistical method for finding the best fits 
is applied. We used a method to compare our model LCs with Fermi LAT data for the Vela 
pulsar: 
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with Vn^i the model (relative flux) value and Y^ \ the measured number of counts (relative 
units) in each phase bin. First we lowered the model LCs resolution, so that both the model 
and data have the same amount of bins N. Next, we smoothed the data using a Gaussian kernel 
density estimator (KDE). The data are treated as being cyclic. For computational efficiency, we 
aligned the maximum peaks of the model and data before calculating x^(a)C)- A contour plot 
of x^ is shown in panel (b) of Figured 


3. Results 
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Figure 2. The sky maps (left) and LCs (right) as predicted from the TPC model using the 
static dipole field for different a and values (deg). 


As an example, we show the sky maps and their corresponding LCs for the TPC model, for 
both the static and retarded dipole fields (Figure [2] and [3|). We used a maximum gap radius of 
Rraa.x = l-2i?LC for both the TPC and OG cases. There are different LCs on the right of each 
sky map corresponding to different C-cuts. In both the figures there appear two dark circles, the 
PCs, followed by two sharp, bright regions near it, called the main caustics, on the sky map. 
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Figure 3. The sky maps (left) and LCs (right) as predicted from the TPC model using the 
retarded dipole field for different a and ( values (deg). 
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Figure 4. Panel (a) indicates our best-fit LC for Vela (see Table [T]). Panel (b) shows the 
contour plot for indicating the best-fit solution. 


The caustic structure is qualitatively different between the two cases, leading to differences in 
the resulting LCs. The caustics seem wider and more pronounced in the retarded dipole case. 
A thin line of emission, due to the ‘notch ’ m is also visible in the latter case. For large a the 
caustics extend over a larger range in for the retarded case compared to the static case. The 
OG models are not visible at all angle combinations and thus do not fill all phase space. This 
is due to emission that occurs below the null charge surface for the TPC model, but not for the 
OG model. The TPC model LCs also exhibit relatively more off-pulse emission. The LCs in 
the OG models are due to emission from only one pole, while both poles are visible in the TPC 
model. 

































































The different model LCs are fitted to the Vela data of Fermi LAT and for each model, we 
constructed a contour plot which indicates the best possible fit. The white marker on the 
contour plot (Figure U panel (b)) indicates the minimum value of log^g and these values 
are shown in Table [TJ The first column shows our different combinations of magnetic field and 
geometric model, the second indicates the minimum value of log^g X^i and the third and fourth 
columns indicate the best-fit a and C from our models. These are for the 5° resolution. We will 
estimate more rigorous errors on these values in future. The fifth column contains the derived (( 
values from the pulsar wind nebula (PWN) torus fitting with a constrained [TS]. The last two 
columns are derived from fits of the rotating vector model to the radio polarization angle (PA) 
versus phase cf) [E]. The best-fit model LC to the Vela LC is shown in Figured panel (a). Our 
best fit is close to values inferred from these independent studies. 

4. Conclusions and future work 

We have studied the effect of different magnetic fields on gamma-ray LC characteristics. We 
utilized the static and retarded vacuum dipole solutions, in combination with the TPC and OG 
geometries. It is evident that the magnetospheric structure and emission geometry determine 
the pulsar visibility and also the 7 -ray pulse shape. We applied our models to the Vela pulsar 
and found a best fit from the OG model using the retarded dipole field, for (a,C) = (80°, 70°). 
This is reasonably close to the value of (a,C) = (75°, 64°) inferred by [16]. In future, we will 
implement an additional magnetic field solution, the offset dipole [ 20 ] and study the effect of 
this solution on the predicted pulsar LCs. 
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